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Abstract 

Electronic transport through a single-molecule magnet Mni2 in a two-terminal set up is calculated 
using the non-equilibrium Green's function method in conjunction with density-functional theory. 
A single-molecule magnet Mni2 is bridged between Au(lll) electrodes via thiol group and alkane 
chains such that its magnetic easy axis is normal to the transport direction. A computed spin- 
polarized transmission coefficient in zero-bias reveals that resonant tunneling near the Fermi level 
occurs through some molecular orbitals of majority spin only. Thus, for low bias voltages, a spin- 
filtering effect such as only one spin component contributing to the conductance, is expected. This 
effect would persist even with inclusion of additional electron correlations. 
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Electronic transport through single-molecule magnets (SMMs) has recently been mea- 
sured in three-terminal setups or using scanning tunneling microscope (STM).— Both 
of the experimental techniques face their own challenges in transport measurements, so do 
theoretical and computational studies of the systems. So far, there were no experimental 
or theoretical inputs on the interfaces between SMMs and electrodes. It is also extremely 
difficult to control these interfaces. SMMs can be differentiated from organic molecules in 
the sense that they have large magnetic moments and large magnetic anisotropy barriers 
caused by spin-orbit coupling. In addition, SMMs differ from magnetic clusters compris- 
ing magnetic elements only because the magnetic ions in SMMs are interacting with each 
other via organic/inorganic ligands. It was shown that the magnetic properties of SMMs 
change significantly with the number of extra electrons added to the molecules,- 1 ^ which may 
greatly impact their transport properties. Several theoretical studies^^ on the transport 
through SMMs, have been, so far, based on many-body model Hamiltonians considering an 
important role of strong correlations in SMMs. However, in these model Hamiltonian stud- 
ies, effects of interfaces and molecular geometries that also play a crucial role in transport 
were not properly included. In this sense, first-principles calculations could complement the 
existing many-body Hamiltonian studies. 
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FIG. 1: (Color online). Extended molecule or scattering region consisting of SMM Mni2 (center) 
attached to Au layers via S atoms and alkane chains. Semi-infinite Au electrodes are considered in 
the calculations (not shown). 

In this work, we present first-principles calculations of transport properties through a 
SMM Mn 12 when the molecule is bridged between Au(lll) electrodes. We identify plausible 
pathways of the transport within the molecule and find a spin-filtering effect in the transport 
caused by the nature of Mn 12 molecular orbitals responsible for resonant tunneling. For 
this study, we use the non-equilibrium Green's function method and the density-functional 
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theory (DFT), with the SIESTA-based 1 ^ quantum transport code, SMEAGOL.— ^ To include 
additional electron correlations, we also add a Hubbard-like U term to our DFT calculations 
using VASP 1 -. Our main results reported in this work do not change with the addition of 
U term except for an increased HOMO-LUMO gap comparable to the experimental value.- 
In our study the following assumptions are made: (1) Despite intermediate transient states, 
the system eventually reaches the stable state. (2) Currents can be obtained from the 
methodology based on the ground-state DFT even when a bias voltage is app lied. (3) 
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we divide 



Interactions of molecules with heat baths are not included. As discussed in Ref. 
the whole system into two parts: (i) bulk left and right Au electrodes, and (ii) a scattering 
region consisting of several Au layers, a SMM Mn 12 , and two linker molecules (S atoms and 
alkane chains). The scattering region [also called extended molecule (EM)] is shown in FigfJJ 
The bulk electrodes are treated semi-infinite and their electronic structures are computed 
using SIESTA. A current is expressed as^ / = (e/h) J dE T(E, V)(f(E - /i L ) - f(E - fi R )), 
where T(E, V) is a transmission coefficient, V is a bias voltage, and f(E — /i L ) and f(E—}j, R ) 
are the Fermi functions for the left and right electrodes with chemical potentials \il and /xr, 
respectively. T(E, V) can be cast as 

T(E, V) = T L (E, V)G\ I (E, V)T R (E, V)G M (E, V) (1) 

where T l R is the density of states of the left or right electrode and Gm is the Green's 
function of the EM. Within the non-equilibrium Green's function method, Gm is solved 
self-consistently in the context of the DFT. 

For spin-polarized DFT calculations with SIESTA we use Perdew-Burke-Ernzerhof (PBE) 
generalized-gradient approximation (GGA)^ for exchange- correlation potential. We gener- 
ate Troullier-Martins pseudopotentials 1 - for Au, Mn, S, O, C, and H with scalar relativistic 
terms and core corrections except for H, as well as basis sets for each element. For Mn 
atoms 3p orbitals are included in the valence states. For an isolated Mni2 the molecular 
orbitals near the Fermi level obtained using the generated basis sets and pseudopotentials 
agree well within 0.1 eV with those using other DFT codes such as NRLM0L 17 > and VASP. We 
also compute the magnetic anisotropy barrier for Mni2 by including the spin-orbit interac- 
tion self-consistently in DFT calculations with a version of SIESTA that includes spin-orbit 
coupling 12 . The barrier we obtain agrees with the experimental valued. All of these tests 
reveal that our generated pseudopotentials and basis sets are good enough to study the 
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FIG. 2: (Color online). Spin-polarized transmission coefficient in zero bias. 



system of interest. 

To construct the bulk leads, we use an optimum bulk lattice constant of 4.166 A, resulting 
in a vertical distance between adjacent layers of 2.405 A. The electronic structure and 
self-energies of Au(lll) electrodes are computed using SIESTA. For the EM we consider a 
geometry in which a Mni2 molecule is oriented such that its magnetic easy axis is normal 
to the transport direction (z axis) and in which six Au layers are included on each side of 
the Mn 12 molecule (Fig. [1]). This geometry completely covers a Mn 12 molecule and prevents 
different Mn 12 molecules from interacting with each other. Mn 12 molecules are not directly 
chemically bonded to Au, and so thiol groups and alkane chains are used to attach Mni 2 to 
Au (Fig. [1]). In the EM, 3 x 3 x 1 k points are sampled. The periodic boundary conditions 
are imposed on the EM along all directions to self-consistently solve for the Green's function 
of the EM. 

The computed spin-polarized transmission coefficient in zero bias is shown as a function 
of energy relative to the Fermi level Ef in Fig. [2j The peaks in the transmission coefficient 
are very narrow due to weak coupling between the electrodes and Mn 12 , as discussed in our 
previous DFT calculation.— Thus, one can identify a one-to-one mapping between individual 
transmission peaks and molecular orbitals. The transmission coefficient from minority-spin 
electrons is zero in the energy region shown in Fig. El The first minority-spin transmission 
peak appears at 0.71 eV above Ef. So contributions to the transmission near Ef are from 
majority-spin electrons only This is because only majority-spin molecular orbitals are near 
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FIG. 3: (Color online). Spin-polarized density of states of the EM projected onto Mn d and C p 
orbitals (majority spin: black, minority spin: red). Minority-spin densities for Mn are zero in the 
energy range shown. The Mn ions are labeled in Fig. 01 The dashed lines (arrows) are for the 
HOMO and LUMO (HOMO-2, HOMO, LUMO+2, LUMO+3) of Mn 12 from the left. The LUMO 
is slightly above the Fermi level. 
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FIG. 4: Positions of Mn ions (filled circles), few O, C, and S atoms. The dashed lines represent 
the alkane chains. The S atoms are directly bonded to the Au leads. 



Ef and they are well separated from minority-spin molecular orbitals. This feature may be 
seen in the spin-polarized projected density of states onto Mn d and C p orbitals (Figs. [3] 
and HJ). The minority-spin densities of Mn d orbitals are zero in the energy range shown 
in Fig. [31 Compared to fourfold symmetric Mn d orbitals in an isolated Mni2, the Mn d 
orbitals in the EM reveal twofold symmetry because of broken symmetry caused by the 
leads (Fig. [3j). By comparing Fig. [3] to Fig. [21 we find that the four transmission peaks 
between E=0 and 0.1 eV are from the majority-spin lowest-unoccupied-molecular-orbital 
(LUMO) and the three levels right above the LUMO (LUMO+n, n=l,2,3), while the four 
peaks between i?=-0.4 and -0.3 eV are from the majority-spin highest-occupied-molecular- 
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orbital (HOMO) and the three levels right below the HOMO (HOMO-n, n=l,2,3). The 
HOMO is mainly from Mn(8) and Mn(12) and the HOMO-1 and HOMO-3 are from Mn(6) 
and Mn(10). The LUMO and LUMO+1 are from Mn(5) and Mn(9) and the LUMO+2 
and LUMO+3 are from Mn(7) and Mn(ll). Although a corresponding transmission peak 
appears at the energy level of each molecular orbital in the range, not all of the orbitals 
would significantly contribute to the current because of the extremely narrow widths of 
peaks. The peaks at _E=-0.335, -0.35, 0.07, and 0.08 eV, have widths of less than 10~ 5 eV, 
and so they would not contribute much to the current. The widths of transmission peaks 
depend on the extent of the electron density overlap or level broadening along the transport 
pathways. For example, for electrons to be transmitted through Mn 12 , they must tunnel 
through the alkane chains. The spin-polarized density of states projected onto the three 
C atoms in each alkane chain shows small peaks at the energies of the HOMO-2, HOMO, 
LUMO+2, LUMO+3 as indicated by the arrows in the bottom panel of Fig. [3j resulting 
in extremely narrow transmission peaks. Thus, more electrons would tunnel through the 
LUMO, LUMO+1, HOMO-1, and HOMO-3. 

For comparison with experimental data, strong correlation effects which are lacking in 
GGA, should be included. Insight on these effects on transport is provided from GGA+U 
studies on an isolated Mn i2 . Our previous GGA+U study on a neutral Mni 2 20 suggests 
that the HOMO-LUMO gap shown in Figs. [2] and [3] would increase so that a higher bias 
voltage would be needed to charge the Mni 2 molecule than appeared in Fig. [2l However, 
with a downward shift of occupied levels and upward shift of unoccupied levels, the main 
features including spin filtering would persist. When the Mni 2 molecule is singly charged 
with inclusion of U, we find that the level corresponding to the LUMO of the neutral Mn 12 
is now filled so that the energy gap for the charged Mni 2 becomes very small. Thus, what 
we learned from the zero-bias case would be still relevant for non-zero bias cases to large 
extent. Further studies on the charged Mni 2 molecule in transport are in progress. 

In summary, we have investigated transport properties through a Mni 2 molecule bridged 
between Au(lll) electrodes using the non-equilibrium Green's function method and spin- 
polarized DFT. We found that a spin-filtering effect occurs in the transmission for non- 
ferromagnetic electrodes. In addition, not all of the molecular orbitals near Ef involved 
with the resonant tunneling, equally contribute to the current. 
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